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ABSTRACT 

We measure the cosmic shear power spectrum on large angular scales by cross-correlating 
the shapes of ~9 million galaxies measured in the optical SDSS survey with the shapes of 
~2.7x 10 5 radio galaxies measured by the overlapping VLA-FIRST survey. Our measurements 
span the multipole range 10 < € < 130, corresponding to angular scales 2° < 6 < 20°. On 
these scales, the shear maps from both surveys suffer from significant systematic effects that 
prohibit a measurement of the shear power spectrum from either survey alone. Conversely, we 
demonstrate that a power spectrum measured by cross-correlating the two surveys is unbiased, 
reducing the impact of the systematics by at least an order of magnitude. 

We measure an if-mode power spectrum from the data that is inconsistent with zero 
signal at the 99% confidence (~ 2.1 cr) level. The odd-parity #-mode signal and the EB cross¬ 
correlation are both found to be consistent with zero (within l<x). These constraints are ob¬ 
tained after a careful error analysis that accounts for uncertainties due to cosmic variance, 
random galaxy shape noise and shape measurement errors, as well as additional errors associ¬ 
ated with the observed large-scale systematic effects in the two surveys. This latter source of 
uncertainty is particularly important for our analysis as it amplifies the errors by a factor ~2.5 
compared to the errors due to cosmic variance, galaxy shape noise and measurement errors 
alone. Our constraints, which probe the power spectrum in the linear regime, are consistent 
with the expected signal in the concordance cosmological model assuming recent estimates 
of the cosmological parameters from the Planck satellite, and literature values for the median 
redshifts of the SDSS and FIRST galaxy populations. 

The cross-power spectrum approach described in this paper represents a powerful tech¬ 
nique for mitigating shear systematics and will be ideal for extracting robust results, with the 
exquisite control of systematics required, from future cosmic shear surveys with the SKA, 
LSST, Euclid and WFIRST-AFTA. 

Key words: gravitational lensing: weak, methods: statistical, cosmological parameters, 
galaxies: distances and redshifts 


1 INTRODUCTION 

Gravitational lensing is the apparent distortion in the shapes of dis¬ 
tant background galaxies due to the path of their light rays being 
altered by the gravitational field of a foreground object. In the weak 
lensing (WL) regime these distortions are so slight that they cannot 
be measured using only one background source; hence WL is an 
intrinsically statistical measurement. WL is considered today to be 
a uniquely powerful tool in astronomy. This is because it can be 
used to estimate the total mass of a foreground lens regardless of 
the constituent matter’s fundamental nature. 

The tiny coherent WL distortions of background galaxy 
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shapes caused by the large-scale structure in the Universe is termed 
cosmic shear. The theoretical basis for cosmic shear was pioneered 
by Gunn (1967), but due to its challenging observational nature 
it was only detected around three decades later (Villumsen 1995; 
Schneider 1998; Bacon et al. 2000; Kaiser et al. 2000; Van Waer- 
beke et al. 2000; Wittman et al. 2000). The field has since flourished 
with numerous experiments delivering results with increasing ac¬ 
curacy (Van Waerbeke et al. 2001; Brown et al. 2003; Bacon et al. 
2003; Fu et al. 2008; Heymans et al. 2012; Jee et al. 2013). Mea¬ 
suring these gravitationally induced distortions at different scales 
provides a wealth of cosmological information. In particular cos¬ 
mic shear studies can constrain the clustering amplitude cr 8 and the 
matter density Q m (Hoekstra et al. 2006; Schrabback et al. 2010; 
Kilbinger et al. 2013). Furthermore with forthcoming and future 
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large surveys containing redshift information (e.g. DES 1 , KIDS 2 , 
HSC 3 LSST 4 , WFIRST-AFTA 5 , Euclid 6 and SKA 7 ), it will be pos- 
sible to measure the evolution of the large-scale matter distribution 
and hence place constraints on the dark energy equation of state w. 

Nearly all cosmic shear studies to date have been performed in 
the optical and near infrared (NIR). The only comparable study in 
the radio band is the analysis of the YLA FIRST 8 survey described 
in Chang et al. (2004). Weak lensing in the radio is still at an early 
stage compared to lensing at optical/NIR wavelengths. The primary 
reason for this is the relatively small number of sources in the radio 
sky at current telescope sensitivities. Current radio surveys can de¬ 
liver source counts of ~ 100 deg -2 , of which ~25% can be used in 
weak lensing studies (Becker et al. 1995; Chang et al. 2004). This 
can be compared to optical surveys which can reach source number 
density of ~10arcmin~ 2 over similar area of sky (e.g. Ahn et al. 
2014). This, together with the inability of radio surveys to deliver 
redshifts for the detected galaxies, means that WL studies in the ra¬ 
dio cannot currently compete with state-of-the-art optical/NIR WL 
experiments. 

This situation is expected to change in the foreseeable fu¬ 
ture. The number counts of radio sources will increase dramati¬ 
cally when radio telescopes breach the l-10juJy sensitivity level, 
at which point radio and optical galaxy counts will be comparable. 
Ongoing and future surveys with the JVLA 9 and e-MERLIN 10 fa¬ 
cilities, the SKA pathfinders MeerKAT 11 and ASKAP 12 , and finally 
the SKA itself, will achieve and ultimately surpass this sensitivity 
level. Additionally with the advent of the SKA, the detection of 
radio galaxies will also be accompanied by estimates of their red- 
shifts from measuring their HI emission. 

Where galaxy numbers are comparable, radio WL studies pro¬ 
vide potential advantages over their optical/NIR counterparts. Ra¬ 
dio interferometers have well-known and deterministic beam pat¬ 
terns. Therefore the instrumental point spread function (PSF) can 
be estimated, at any position on the sky, to very high accuracy. This 
is in contrast to optical studies where the PSF is typically mapped 
across the sky using point sources (stars) and is then extrapolated 
to the positions of the galaxies. This PSF measurement and inter¬ 
polation scheme will inevitably introduce errors to the subsequent 
galaxy shape measurements at some level. 

Future radio surveys are also expected to be sensitive to the 
star-forming galaxy population to higher redshifts than is typi¬ 
cally probed with optical surveys (Brown et al. 2013, 2015). Ra¬ 
dio WL surveys are therefore expected to probe the Universe at 
earlier times. Additionally, Brown & Battye (2011) demonstrated 
the potential use of polarisation information (which comes rela¬ 
tively easy with radio observations) to estimate the orientations of 
galaxies prior to lensing. This technique can be used to suppress the 
impact of intrinsic galaxy alignments, a key astrophysical system¬ 
atic effect that contaminates weak lensing studies at all wavelenghts 


1 Dark Energy Survey, http: //www. darkenergysurvey. org 

2 Kilo Degree Survey, http: //kids. strw. leidenuniv. nl 

3 Hyper Suprime-Cam, http: //www. naoj . org/Proj ects/HSC 

4 Large Synoptic Survey Telescope, http: //www .lsst.org 

5 Wide-Field IR Survey Telescope, http: //wfirst. gsfc. nasa. gov 

6 Euclid satellite, http: //sci . esa. int/euclid 

7 Square Kilometre Array, http: //www. skatelescope. org 

8 Very Large Array; Faint Images of the Radio Sky at Twenty centimetres 

9 Jansky VLA, https://science.nrao.edu/facilities/vla 

10 http://www.e-merlin.ac.uk 

11 http://www.ska.ac.za/meerkat/ 

12 http://www.atnf.csiro.au/proj ects/askap/index.html 


(see Joachimi et al. 2015, Kiessling et al. 2015 and Kirk et al. 2015 
for a recent overview). Radio polarization measurements can po¬ 
tentially also be used to significantly reduce the number of galaxies 
required for a weak lensing study. Looking to even higher redshifts, 
Pourtsidou & Metcalf (2014) have shown that one can exploit 21cm 
SKA surveys to perform weak lensing studies without having to re¬ 
solve or even identify individual galaxies. For more details on the 
SKA WL prospects in general see Brown et al. (2015). 

Another key advantage, which we aim to exploit in this paper, 
arises when radio and optical/NIR data are combined. WL anal¬ 
yses rely on extremely accurate (and, to a lesser extent, precise) 
measurements of galaxy shapes. However instrumental systematic 
effects can introduce correlated errors in the measurements. Since 
it is primarily correlations in the galaxy shapes that one wishes 
to measure for cosmology, such instrumental systematics repre¬ 
sent a significant problem for WL surveys. The source of sys¬ 
tematic effects in optical/NIR telescopes and radio interferome¬ 
ters are of a very different origin and nature, and are therefore 
not expected to correlate. By cross-correlating the shape informa¬ 
tion from two such experiments, one therefore expects to cancel 
out any instrumental systematics and obtain an unbiased cosmic 
shear measurement (Jarvis & Jain 2008; Patel et al. 2010). In this 
paper, we aim to demonstrate the potential benefits of such cross- 
correlaton techniques by performing an exploratory optical-radio 
cross-correlation weak lensing analysis of the VLA FIRST survey 
(Becker et al. 1995) and the Sloan Digital Sky Survey (SDSS, York 
et al. 2000). These two surveys are well suited for this investiga¬ 
tion as they provide approximately 10,000 deg 2 of overlapping sky 
coverage. In contrast to most other cosmic shear analyses, our mea¬ 
surements probe the shear power spectrum on large scales, at multi¬ 
poles £ < 100, corresponding to angular scales above a few degrees. 
We perform our analysis using a cross-power spectrum approach, 
which provides for a natural decomposition of the signal into its 
constituent E- and 5-modes, and also fully accounts for the sur¬ 
veys’ masks and geometries on the curved sky. 

The paper is organised as follows. In Section 2 we present 
the relevant weak lensing formalism, and the details of our algo¬ 
rithm for estimating the cosmic shear cross-power spectrum from 
the data. Section 3 describes the two surveys that we use in the 
analysis. In Section 4 we describe how we construct shear maps 
from the two datasets. We identify large-scale systematic effects in 
these maps which dominate the measured shear auto-correlations 
in either dataset alone. We demonstrate using simulations in Sec¬ 
tion 5 that the systematics can be removed by cross-correlating the 
shape information from the two datasets. In Section 6 we present 
the results from the analysis of the real data where we measure the 
shear cross-power spectrum from the FIRST and SDSS datasets. 
We discuss our results and conclude in Section 7. 


2 WEAK LENSING THEORY 


Distortions in the shapes of distant galaxies can be quantified via 
measurements of their ellipticity, which is a spin-2 quantity: 



The two components of this matrix represent the two orthogonal 
modes of an ellipse: e\ describes elongation or compression along 
or perpendicular to a user-defined reference axis while e 2 describes 
elongation or compression at an angle n/4 to the reference axis. 
Galaxy ellipticity can also be described in terms of a semi-major 
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axis (a), a semi-minor axis ( b ), and a position angle ( 0 ) as 

6\ = £ x cos( 20 ), ( 2 ) 

£2 = £ x sin(20). (3) 


Here, e = ( a-b)/(a+b ) is the modulus of the galaxy’s ellipticity and 
0 is the position angle of the galaxy with respect to the reference 
axis of the chosen coordinate system. 

In the absence of intrinsic galaxy alignments, galaxy ellipticity 
is an unbiased estimator for the lensing shear field, which is also a 
spin -2 quantity, 


<e, 7 > = 7ij = 


n 

72 



(4) 


where the angled brackets denote an ensemble average. Since grav¬ 
ity is a potential theory, the shear at angular position Q can be re¬ 
lated to a lensing potential ( 0 ) as 

yij(£l) = \6i6j - I(5£(5 2 j^(n), (5) 

where = r(Sij - r^V;) is a dimensionless, transverse differential 
operator, and 5 2 = 6i6 j is the transverse Laplacian. The lensing 
potential can in turn be related to the 3-D gravitational potential, 
®(r) by (e.g. Kaiser 1998) 



where r is the comoving distance to the sources. 

Note that the ellipticity and shear fields of equations (l)-(4) 
are defined with respect to an arbitrarily defined reference axis. In 
the case where one is interested in the weak lensing distortion in a 
population of background sources due to the presence of a known 
foreground object (or “lens”), it is more natural to consider the tan¬ 
gential and rotated shear (or ellipticity), defined by 


7t = 71 cos(20) + y 2 sin(20), (7) 

y r = -yi sin( 20 ) + y 2 cos( 20 ), ( 8 ) 


where 6 is the position angle formed by moving counter clockwise 
from the reference axis to the great circle connecting each source- 
lens pair. The tangential shear, y f , describes distortions in a tan¬ 
gential and/or radial direction with respect to the lens position. The 
rotated shear, y r , describes distortions in the orthogonal direction, 
at an angle ±/r/4 to the vector pointing to the lens position. 

As indicated by equation (4), the shear components can be 
estimated by averaging over a set of observed background galaxy 
ellipticities. The shear field can itself be related to the lensing con¬ 
vergence field (which measures the projected surface mass density 
contrast) by (Kaiser & Squires 1993) 

k = d~ 2 didj7ij, (9) 

where d ~ 2 is the inverse 2-D Laplacian operator defined by 

<T 2 = — f d 2 r' In \r - r'\. ( 10 ) 

2 n J 

A second decomposition of the shear yields the odd-parity diver¬ 
gence field, 

p = d- 2 e?‘ i d j d n y ij , ( 11 ) 


where e? is the Levi-Civita symbol in two dimensions, 



( 12 ) 


Since gravitational lensing produces no /?-modes, at least to first 
order, the divergence field is a useful quantity for testing systematic 
effects in weak lensing studies. 


2.1 Weak lensing in spherical harmonic space 

As mentioned above, lensing shear is a spin-2 field, i.e. it trans¬ 
forms as 7 —> ye l2(f> under rotation by 0 (where, for convenience 
we have defined the complex shear y - y x + iy 2 ). One may thus 
expand y and its complex conjugate in terms of the spin-weighted 
spherical harmonics, s Y {m (Newman & Penrose 1966) as 


y(tl) = y,(Q) + iy 2 (ti) 


~ "*■ ) ]Y I'm (^2), 

(13) 

£m 


r*(fl) = ym-iy 2 hh 


= YY Kem ~ -2 

(14) 


£m 


where 5 denotes the spin and the summation in m is over -l < 
m < l. K£ m and Pt m are the spin -2 harmonic modes of the so-called 
electric (i.e. “gradient” or “£-mode”) and magnetic (i.e. “curl” or 
“5-mode”) components of the shear field respectively, which we 
further identify as the harmonic space versions of the lensing con¬ 
vergence and divergence fields mentioned in the previous section. 
Using the orthogonality of the spin-s spherical harmonics over the 
whole sphere, 

J*d£l s Yt m (Q) s Yp m , (Q) = 6a f d mm >, (15) 

where the spin states must be equal, the harmonic modes of the k 
and p fields can be found directly from the shear field, 7 ; 

Ktm = \fd£l[ 7 (£l) 2 r (m (£l) + y(£l)- 2 Y* tm m, (16) 

p lm = y Jdti [ 7 (n) 2 r;,(n)- y *(a)_ 2 ^,(n)]. (17) 

Taking the average values over the sphere of products of the har¬ 
monic coefficients of the k and p fields, one can construct three 
possible power spectra: 

Cp - + l yi Kim K Zm' (18) 

m 

~ 21 + 1 PtmPtm' (19) 

m 

C ?= 2 €TT (20) 

m 

The parity invariance of weak lensing suggests that and Cf 
should be equal to zero in the absence of systematics. However, 
systematic effects (both instrumental and astrophysical) can give 
rise to a non-zero . Finite fields and boundary effects can also 
lead to leakage of power between the three spectra. 

In the limit of weak lensing, the two-point statistical properties 
of the shear and convergence fields are the same (Blandford et al. 
1991) so that Cy = Cf. Finally, we can relate the convergence 
power spectrum to the 3-D matter power spectrum P§(k, r ) through 
(e.g. Bartelmann & Schneider 2001) 
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where we have assumed a flat Universe, and we have now con¬ 
sidered the more general case where one measures the cross¬ 
correlation signal between two different surveys, here denoted with 
indices i and j. In equation (21), H 0 is the Hubble constant, a is 
the scale factor of the Universe, r is comoving distance, r H is the 
comoving distance to the horizon and Q m is the matter density. The 
weighting, W t (r), is given in terms of the normalised source distri¬ 
bution for each survey, Gi(r)dr = pi(z)dz: 

( 22 ) 

2.2 Power spectrum estimation on a cut sky 

The power spectrum estimators of equations (18)-(20) are only un¬ 
biased in the case of full sky coverage and in the absence of noise. 
However, in the analysis which follows we shall be estimating these 
power spectra from noisy data covering only a fraction of the sky. In 
order to correct for the effects of limited sky coverage and noise, we 
adopt a “pseudo-Q” analysis, originally developed in the context 
of CMB polarization observations. We provide a brief summary of 
the technique here and refer the reader to Brown et al. (2005) for 
further details. 

In the presence of finite sky coverage and/or a mask that ex¬ 
cludes certain regions of the sky, one can estimate the so-called 
pseudo-Q’s directly from a set of shear maps (using fast spheri¬ 
cal harmonic transforms from e.g. the HEALPIX software 13 , Gorski 
et al. 2005) as: 

G (W) = X (23) 

m 

where {A, B } are any of the three possible combinations of k and p 
and {/, j] denote the two different surveys as before. In the following 
sections, we will mostly be interested in taking the cross-power 
spectra of the FIRST and SDSS datasets, i.e. where i denotes the 
FIRST dataset and j denotes SDSS. When cross-correlating two 
datasets, there are then in principle two possible Cf power spectra. 
For simplicity, we will take the average of these and present a single 
Cf cross-power spectrum. In equation (23), A im and B im denote the 
spin-2 spherical harmonic modes (found using equations 16 and 17) 
of weighted versions of the lensing shear field, 

r(D) = W(Q)y(Q), (24) 


Wi(r) : 


f 


dr' G/(re¬ 


calculate the Mw matrix exactly from the pseudo-Q cross-power 
spectra (estimated using equation 23) of the weighting functions 
used for the two different surveys, W;(Q) and Wj(Cl). Further details 
and explicit formulae for calculating M u > are provided in Brown 
et al. (2005). 

Recovering an estimate of the power spectra on the full sky 
is then only a matter of inverting equation (25). When calculating 
auto-power spectra (from a single survey), one also needs to include 
a correction for the noise bias arising from the intrinsic dispersion 
in galaxy shapes and measurement errors. For practical reasons, it 
is often also convenient to estimate the power spectrum in terms of 
“band powers” which recover the average power across a range of 
multipoles, AG Including these features, the power spectrum esti¬ 
mator becomes 

= X K '*’ X Ob f e(Ct - (AQmcX (26) 

b' l 

where P b denotes the debiased combined estimator for bandpow- 
ers of the Cf, Cf and Cf power spectra. In equation (26), Q are 
the pseudo-Q’s measured from the data and (AQ are the noise-only 
pseudo-Q’s estimated from a suite of Monte Carlo (MC) simula¬ 
tions which contain realisations of the noise, and also incorporate 
the effects of the survey masks. Note that, for our cross-power spec¬ 
trum analysis, the noise is_uncorrelated and so the noise subtraction 
step (subtraction of the (AQmc term) is not required, although we 
have nevertheless implemented it in our analysis pipeline for com¬ 
pleteness. 

K b b f is the binned coupling matrix, constructed from Mw as 

Kbb' = ^ O b £ ^ Mw F > Qi ’ b >, (27) 

t v 


where O b i is a binning operator that bins the Q’s into band powers 
and Qib is the inverse operator which “unfolds” the band powers 
into individual Q’s. The function, Q describes the smoothing ef¬ 
fect on the underlying shear fields due to the fact that the sky maps 
are pixelized. This function is provided by the healpix software 
that we use to pixelize the sky and perform the spin spherical har¬ 
monic transforms. 

The binning operator we choose to use is 


Obi - 


1 i(i+ 1) 

2 n/b+l)_Jb) ’ 

1 low 1 low 

o, 


otherwise , 


(28) 


where W(Q) is an arbitrary weighting function which can be used 
to exclude certain regions of the sky and/or downweight noisier sky 
pixels. 

Grouping all three power spectra into a single vector, Q = 
{C f, Cf, Cf }, one can show that the pseudo-Q’s measured from 
the weighted shear fields are related to the true power spectra on 
the full sky via 

C ( = Y J M U ’Ce, (25) 

r 

where we have introduced the coupling matrix, M U r. This cou¬ 
pling matrix fully encodes how the survey geometries and masks 
(described via the function, W(D)) mix modes, both within a sin¬ 
gle spectrum (e.g. Cf —> Cf), and also between spectra (e.g. 
Cf —> Cf), and it can be calculated exactly for any given set of 
weighting functions, W(f2). Specifically, in our case where we are 
estimating the cross-power spectra between two datasets, we can 

13 see http://healpix.sourceforge.net 


where iff denotes the lower edge of band b. With this choice, the 
quantity i(i+l)C^/2n is approximated as constant within each band 
power. 

3 THE SURVEYS 
3.1 FIRST data 

The VLA FIRST survey (Becker et al. 1995) is a 10,575 deg 2 sur¬ 
vey of the sky in the declination range -10° < 6 < 70°, and right 
ascension range 8hrs < RA < 17hrs. Of the survey’s total sky 
coverage, 8,444 deg 2 are in the north galactic cap and 2,131 deg 2 
are in the south. The data were gathered using the VLA at L-band 
(1.4 GHz), in B-configuration (maximum baseline length of 11.1 
km) and in snapshot mode (<5 minutes per observation). The ob¬ 
servations were conducted between 1993 and 2011. The coverage 
for the southern cap is discontinuous due to poor weather and sys¬ 
tem failures during the 2011 observations. 
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The beam of the telescope at 6 > +4°33' can be approxi¬ 
mated by a circular Gaussian with a full width at half maximum 
(FWHM) of 5.4”. Below 6 - +4°33' the beam is elliptical with 
a FWHM of 6.4”x5.4”, with the major axis running north-south. 
For 21hrs < RA < 3hrs and 6 < -2°30' the ellipticity of the 
beam increases further to 6.8”x5.4”. The survey achieved an RMS 
noise level of -0.15 mJy, and delivered positions, flux densities and 
shape information for -1 million radio sources, of which -40% are 
resolved. Shape information was extracted by fitting an elliptical 
Gaussian model to each source. The FIRST catalogue provides this 
shape information via a major and minor axis and a position angle 
(PA) for each object. Note that the FIRST analysis pipeline includes 
a deconvolution step to remove blurring due to the telescope’s PSF. 
More information regarding the contents of the FIRST catalogue 
and how it was generated is available from the FIRST website 14 . 

Mainly due to the relatively sparse coverage in visibility (or 
“wv”) space - an unavoidable result of the observation strategy - 
the deconvolution/CLEANing 15 of the FIRST sources was imperfect. 
The catalogue therefore contains information about the possibility 
of a detected source being a residual sidelobe of a nearby brighter 
source. Of the total number of FIRST sources, 76.4% had a proba¬ 
bility of being spurious, P(S), less than 5%. For 17% of the sources 
the probability was between 5 and 50%. The remaining 6.5% had 
more than a 50% chance of being a sidelobe of a nearby bright 
galaxy. Additional information about the sources was generated by 
cross-correlating their positions with positions of optical sources 
from the SDSS and 2-micron All-Sky Survey (2MASS, Skrutskie 
et al. 2006) catalogues. 

The FIRST data have previously been used to make a detec¬ 
tion of cosmic shear in the radio band by extracting source shape 
information directly from the uv plane (Chang et al. 2004). Those 
authors measured the 2-point correlation function (2PCF) on an¬ 
gular scales, 1° < 6 < 40°. On scales 1° < 6 < 4° the measured 
/5-modes where consistent with zero while a significant signal was 
detected in the k- modes at the -3 cr level. From their results, Chang 
et al. (2004) placed a joint constraint on the clustering amplitude, 
cr 8 and the median redshift of the FIRST sources without an op¬ 
tical counterpart, z m of cr 8 (z m /2 ) 0 - 6 ^ 1.0 + 0.2. Adopting a prior 
of cr 8 = 0.9 + 0.1 they obtained the constraint z m - 2.2 + 0.9 
( 68 % CL). This value is consistent with existing models of the ra¬ 
dio source luminosity function. The results are encouraging as they 
show that although radio surveys lack the sensitivity at the moment 
to compete with optical lensing surveys, they can deliver shape in¬ 
formation which, with the proper processing, can be used for weak 
lensing studies. 

For the purposes of this study we acquired the 13JUN05 version 
of the FIRST catalogue. The Chang et al. (2004) study estimated 
the median redshift of the entire FIRST sample to be between 0.9 
and 1.4 (depending on the model of the radio source luminosity 
function used). Using the SKA simulated skies (S 3 ,Wilman et al. 
2008) - a computer simulation based on more recent data, and de¬ 
signed to model the radio and submillimetre Universe - we estimate 
the median redshift of the entire FIRST sample to be z m -1.2. 


14 FIRST website: http: //sundog. stsci . edu/index. html 

15 The CLEAN algorithm (Hogbom 1974) is a standard technique com¬ 
monly used in radio astronomy to deconvolve images for the effects of a 
finite PSF. The algorithm models the data as a collection of point sources 
and, starting with the brightest source, iteratively subtracts (from the uv- 
data) the flux associated with each source in order to detect and characterise 
fainter objects. 


3.2 SDSS data 

The SDSS is an ongoing (since 2000) optical survey of the north 
and south galactic caps north of declination -15°, covering - 14,500 
deg 2 of the sky. The survey uses a 2.5 metre telescope located at 
Apache Point Observatory (APO) in the Sacramento mountains in 
south New Mexico (Abazajian et al. 2003; Ahn et al. 2014). In 2008 
the experiment entered a new phase called SDSS-III in which new 
instruments came into operation (Eisenstein et al. 2011). One of 
the SDSS-III surveys, the Baryon Oscillation Spectroscopic Sur¬ 
vey (BOSS; Dawson et al. 2013), obtained spectra of -1.4 million 
luminous red galaxies (LRGs) and -0.3 million quasars which have 
subsequently been used to place stringent constraints on cosmolog¬ 
ical models (Aubourg et al. 2014). 

The SDSS data have been used in numerous weak lensing ex¬ 
periments probing different angular scales on the sky. The data have 
proved particularly useful for galaxy-galaxy lensing (e.g. Fischer 
et al. 2000; Mandelbaum et al. 2006) and cluster lensing (e.g. Man- 
delbaum et al. 2008; Sheldon et al. 2009; Rozo et al. 2010). In 
terms of cosmic shear measurements, Lin et al. (2012) and Huff 
et al. (2014) have both exploited the deep observations of SDSS 
“Stripe 82” to perform cosmic shear analyses. Both of these latter 
studies detected a significant /c-mode signal on degree scales with 
/5-modes found to be consistent with zero. Derived constraints on 
Q m and cr 8 were also found to be consistent between these two in¬ 
dependent studies. 

In this study, we have used the tenth data release (DR 10) of 
SDSS (Ahn et al. 2014). DR 10 delivers photometric information 
across 5 bands (u, g, r, i, z) for -500 million galaxies and stars (the 
PHOTO sample), and spectra for -2.5 million of them (hereafter, 
the SPECTRO sample). The mean magnitude across all PHOTO 
galaxies and bands is -22.5. The median magnitude and redshift for 
the SPECTRO galaxy population is -21 and ~0.3 respectively. The 
SPECTRO sample includes two distinct populations one of which 
peaks at z - 0.15 and the second at z - 0.5. Only a few galaxies in 
the sample are at redshifts greater than - 0.7. We chose to download 
all sources from the SDSS DR 10 PHOTO sample that were identi¬ 
fied in i-band as galaxies with magnitudes i mag < 26. For the -38.5 
million sources that met the criteria we acquired information about 
their positions (RA and 6) and multi-band information describing 
their flux densities, galaxy types, second and fourth moments, the 
reconstructed PSF at each galaxy position, the galaxy magnitudes 
and quoted errors, and finally their extinction. 


4 SHEAR MAPS AND TESTS FOR SYSTEMATICS 

Telescope induced effects can severely bias a WL study. Therefore 
systematics must be accounted for and if possible corrected. In this 
section, we describe how we construct pixelized maps of the FIRST 
and SDSS shear fields from the two galaxy catalogues, and we as¬ 
sess the resulting maps for systematic effects. 

4.1 FIRST analysis 

4.1.1 Source selection 

A detailed study of the systematic effects in the FIRST survey was 
performed in Chang et al. (2004). To control systematics, Chang 
et al. (2004) first performed a series of cuts on the source catalogue 
in order to remove from the sample those sources whose shapes 
were most likely to be corrupted due to residual systematic effects. 
In this analysis, we have applied as many of their source selection 
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this test clearly indicate the presence of residual systematics in the 
FIRST galaxy ellipticities, which are strongly correlated with the 
known PSF of the VLA-FIRST observations. These residuals are 
clearly spurious, and are almost certainly the result of an imper¬ 
fect deconvolution and/or cleaning of the FIRST data during the 
imaging step in the FIRST data reduction. In the analysis which fol¬ 
lows, we will attempt to include the effects of these residual elliptic- 
ity correlations by incorporating a model of the systematic in our 
MC simulations. However, we also expect our cross-power spec¬ 
trum approach - whereby we cross-correlate the FIRST and SDSS 
galaxy shapes - to be robust to the presence of this FIRST-specific 
systematic effect. 


4.1.3 FIRST shear maps 


Figure 1 . Results of the test for residual beam systematics in the FIRST 
galaxy shapes, as described in Section 4.1.2. A large radial (negative y t , 
black triangles) distortion, decreasing as a function of angular scale, 6, is 
found in the FIRST galaxy ellipticities when their shapes are stacked around 
the FIRST galaxy positions. The rotated shear (y r , red circles) is consistent 
with zero. 


criteria as possible given the information that was available to us in 
the FIRST catalogue. We therefore discarded sources that had an 
unresolved deconvolved minor axis, or a deconvolved major axis 
that was greater than 7 arcsec or smaller than 2 arcsec. We also re¬ 
moved sources that had an integrated flux that was smaller than 1 
mJy. Finally we have included only sources that had a possibility 
of being a sidelobe of a nearby brighter source P(S) <5%. After 
applying these cuts we were left with ~2.7xl0 5 radio sources (from 
an initial sample of ~1 million catalogue entries). For all remain¬ 
ing sources, we formed the e\ and e 2 ellipticity estimates from the 
catalogued major and minor axes, and position angle, according to 
equations (2) and (3). 


4.1.2 Residual systematics in FIRST ellipticities 

We have performed a number of quality checks on the data immedi¬ 
ately after applying the cuts described above. By far the most infor¬ 
mative of these was a test for corruption of the galaxy ellipticities 
due to residual sidelobe contamination resulting from an imperfect 
deconvolution and/or cleaning procedure. To perform this test, we 
have made use of the g-g lensing tangential and rotated shear con¬ 
structions ( y t and y r ; equations 7 & 8) and have stacked the shapes 
of the selected FIRST galaxies around the positions of all of the 
sources in the original FIRST catalogue (including both resolved 
and unresolved objects). The results of this test are shown in Fig. 1 
where we see a strong distortion in the FIRST galaxy shapes ori¬ 
ented radially from the central stacking position. This signal per¬ 
sists when we randomly choose a sub-set of the FIRST sources’ 
shapes and/or positions. Moreover, the signal shows no obvious de¬ 
pendency on the flux of the stacked and/or central sources. 

To further investigate the origin of this signal we repeated the 
stacking analysis in (ARA, AS) space (i.e. where the positions of all 
of the central sources are re-mapped to ARA = 0; A 6 = 0). The 
resulting maps of tangential and rotated shear are shown in Fig. 2. 
The 6-arm star pattern apparent in the y t map, displayed in the left 
panel of Fig. 2, closely matches the synthesised beam (or PSF) of 
the VLA in “snapshot mode” which was the observation mode em¬ 
ployed during collection of the FIRST survey data. The results of 


The pseudo-Q power spectrum analysis outlined in Section 2.2 re¬ 
quires as input pixelized maps of the shear components. We create 
these maps by simple averaging of the ellipticities of the FIRST 
sources which fall within each sky pixel and we use the healpix 
scheme (Gorski et al. 2005) to define the sky pixelisation. We 
choose the pixel size according to the density of sources in the 
FIRST survey: with ~ 2.7x10 5 sources surviving the cuts described 
in Section 4.1.1, a pixel size of ~1 deg 2 results in a mean pixel 
occupation number of ~20 galaxies. This choice ensures that our 
maps are contiguous over the survey area and are mostly free from 
unoccupied pixels within the survey boundaries. We therefore set 
the healpix resolution parameter, - 64, which corresponds to 
a pixel size of ~0.85 deg 2 . 

Fig. 3 shows maps of the shear components, y\ and y 2 , as re¬ 
constructed from the FIRST galaxy ellipticities. These maps are 
dominated by random noise due to the intrinsic dispersion in galaxy 
shapes and measurement errors. Visual inspection of the figure does 
also suggest a weak dependence of the y\ amplitudes on declination 
though some of this may simply reflect the varying galaxy number 
density, which also shows some dependence on declination. (Later 
in Section 6, we perform a suite of null-tests on our cross-power 
spectrum results, one of which is to split the datasets according 
to declination. This test should highlight any significant problems 
associated with any declination-dependent systematic effects that 
persist in the FIRST shear maps.) 

In addition to the shear maps, we construct a weight map (the 
W(Q) of equation 24) as simply the number of galaxies that lie 
within each pixel. In the case of uniform measurement errors on the 
ellipticities, setting VF(D) to be equal to the galaxy number density 
is equivalent to an inverse-variance weighting scheme, which, for 
the noise-dominated regime in which we are working, is the opti¬ 
mal weighting scheme to use. To limit our sensitivity to individual 
outliers in the galaxy ellipticity distribution, we set W(f2) = 0 for 
pixels that contain less than 5 galaxies. Fig. 4 shows the weight 
map that results from this process. We use this map to weight the 
FIRST shear data in the subsequent power spectrum analysis. 


4.2 SDSS analysis 

4.2.1 Shape measurements and initial source selection 


The SDSS catalogue provides shape information for each detected 
source as ellipticity components derived from second moments, 


Qxx Qyy 
Qxx F Qyy 


(29) 
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Figure 2. Maps of the tangential shear (y f , left panel ) and rotated shear (y r , right panel) as a function of the separation in RA and 8 from the central stacking 
positions. The tangential shear map reveals a negative amplitude 6 -arm star pattern which closely resembles the structure in the synthesised beam (PSF) of the 
VLA snapshot observations. 



Figure 3. Maps of the y\ (left panel) and 72 (right panel ) shear components, constructed by simple averaging of the FIRST galaxy ellipticities within each 
pixel. 



Figure 4. The galaxy number density map used to weight the FIRST shear 
field (Fig. 3) in the power spectrum analysis. 


2gxy 

Qx x “f Qyy 


(30) 


where the adaptive moments are measured directly from the SDSS 
images according to 


Qtj = 


i 


d£lI(Q)w(Q)QiQj. 


(31) 


Here, I(Q ) is the object intensity at position Q = (Q x , Q y ) and {/, jj 
can take values {x, x }, {y,y} or {x,y}. Ideally, adaptive moments are 
measured using a radial weight function, w(H), which is iteratively 
adapted to the object’s size. In practice, the SDSS pipeline used 
a Gaussian weighting function with a size matched to that of the 
object being fitted. Such an approach has been shown to be nearly 
optimal (Bernstein & Jarvis 2002) and is much quicker. The SDSS 
catalogue’s galaxy ellipticities need to be corrected for PSF smear¬ 
ing and atmospheric seeing. To facilitate this, the catalogue also 
provides moment information that describes the PSF reconstructed 
at each galaxy’s position. We follow the procedure advocated by the 
SDSS team which corrects the catalogue ellipticities for the effects 
of seeing and PSF anisotropy according to 

e cor r = (f meas _ _ R ) 

e corr = {e meas _ Re P*fy( l _ (32) 

where e™ ms are the galaxy ellipticity components, e p F are the PSF 
ellipticity components reconstructed at each galaxy position and R 
is a parameter formed from a combination of the second and fourth- 
order moments of the galaxy and PSF shapes. (For further details, 
and an explicit formula for R, see the description on the SDSS 
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Figure 5. Maps of the yi {left panel ) and 72 (n'g/zt panel) shear components, constructed by simple averaging of the SDSS galaxy ellipticities within each 
pixel. These maps were constructed from the ~25 million galaxies remaining in the SDSS shear catalogue immediately after the PSF correction and source 
selection steps described in Section 4.2.1 were implemented. Significant large-scale systematic effects are clearly evident in these maps. 


DR 10 catalogue webpage 16 .) As mentioned in the SDSS catalogue 
description, the correction of equation (32) is not exact and results 
in a small bias (Bernstein & Jarvis 2002; Hirata & Seljak 2003). 
However, rather than implementing a more sophisticated correction 
algorithm, we will mostly appeal to the cross-correlation approach 
of Section 2 to limit our sensitivity to any residual biases that re¬ 
main in the shear catalogue after source selection and implement¬ 
ing the PSF correction of equation (32) (although see below for an 
additional source selection cut that we have applied based on the 
strength of the reconstructed PSF anisotropy, prior to construction 
of our final shear maps). 

To limit the impact of very poorly measured and/or noisy 
shape estimates, we perform an initial source selection on the SDSS 
galaxy sample. To select galaxies, we follow the works of Bern¬ 
stein & Jarvis (2002) and Mandelbaum et al. (2008). We thus ex¬ 
clude sources for which the modulus of the corrected ellipticity 
6 = ( 6 j 2 + > 4 or if the uncertainty in the corrected modulus, 

cr £ > 0.4. We also only retain sources that have an i-band extinc¬ 
tion that is less than 0.2. In addition, we retain only sources that 
have both a PSF and DeVaucouleur r-band magnitude greater than 
22, and both a PSF and DeVaucouleur i-band magnitude greater 
than 21.6. Finally, to exclude very small galaxies, we only include 
sources with Rf > 1/3 where the resolution factor, Rj is given by, 


R f = l- 


QPSf + QPSf 
\cxx ' is iyy 

Qxx T Qyy 


(33) 


With these selection criteria applied, ~25 million SDSS sources 
remain in the catalogue. 


4.2.2 Shear systematics and additional source selection 

Using the shape information of the remaining SDSS sources we 
construct maps of the shear components 71 and 72 on a V side = 64 
healpix grid using simple averaging, the same as was used in the 
previous section for the FIRST maps. The resulting maps are shown 
in Fig. 5. Note that, because of the large increase in galaxy num¬ 
bers, the noise in much lower in the SDSS maps than in the FIRST 
maps of Fig. 3. Visual examination of Fig. 5 shows that these initial 
SDSS shear maps are clearly contaminated with large-scale spuri¬ 
ous arcs. Comparing the structure and morphology of these large- 
scale systematics with a map of the reconstructed PSF anisotropy 


16 http://www.sdss3.org/drl0/algorithms/classify.php 


reveals a strong correlation, indicating that the observed systemat¬ 
ics are residual biases that the approximate correction scheme of 
equation (32) has failed to account for. 

In order to limit the impact of this residual bias, we have im¬ 
plemented an additional cut on the galaxy catalogue based on the 
strength of the PSF anisotropy at each galaxy position. We discard 
all remaining SDSS sources for which the reconstructed PSF ellip¬ 
ticity modulus is e > 0.08. (We have experimented with this thresh¬ 
old value and have found 0.08 to be the minimum value which still 
suppresses the systematics to approximately the level of the ran¬ 
dom noise.) Implementing this additional cut, we were left with ~9 
million SDSS sources. The SDSS shear maps constructed from this 
revised catalogue are shown in Fig. 6 . Visual inspection of the new 
maps demonstrate a large reduction in the amplitude of large-scale 
systematics, which now appear to be at or below the level of the 
random noise. 

This additional cut on the galaxy catalogue increased the noise 
in the SDSS shear maps by a factor ~ V3. However, the measured 
SDSS auto shear power spectrum and the uncertainties associated 
with the large-scale systematics were reduced by a factor of ~ 10 . 
Although the FIRST-SDSS cross-power spectrum analysis will be 
robust to this SDSS-specific systematic effect, suppressing the PSF 
systematic to the level of the noise at this stage in the analysis 
is still beneficial in terms of the error performance of the cross¬ 
power spectrum analysis: without this additional selection based on 
the PSF anisotropy amplitude, the uncertainties in the cross-power 
spectrum increase by a factor of ~ 5 (see Section 6 for more details). 

Finally, we generate a weight map for the SDSS survey using 
the same method as was used in the previous section for the FIRST 
data. The resulting weight map is shown in Fig. 7 and shows a non- 
uniform distribution of SDSS source number density. In particular, 
the very deep coverage of Stripe 82 is clearly visible as an excess 
of galaxies near the Galactic anti-centre, at declinations 8 ~ 0°. 


5 SIMULATIONS 

In this section, we describe our approach to creating simulations of 
the SDSS and FIRST shear datasets. Our simulations serve several 
purposes. Firstly, the pseudo-Q power spectrum estimator (equa¬ 
tion 26) includes a noise-subtraction step which debiases a mea¬ 
surement of the shear power spectrum for the effects of random 
shape noise (due to both intrinsic galaxy shape noise and measure¬ 
ment errors). For this purpose, we estimate the noise bias, (A^)mc 
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Figure 6. Maps of the yi (left panel ) and 72 (right panel ) shear components, constructed by simple averaging of the SDSS galaxy ellipticities within each 
pixel. These maps were constructed from the ~9 million galaxies remaining in the SDSS shear catalogue immediately after the additional source selection 
based on the strength of the PSF, described in Section 4.2.2, was implemented. The large-scale systematic effects apparent in Fig. 5 have been substantially 
reduced with this additional cut on the galaxy catalogue. 



Figure 7. The galaxy number density map used to weight the SDSS shear 
field (Fig. 6 ) in the power spectrum analysis. 


as the average pseudo-Q power spectra measured from a suite of 
simulations containing realisations of the random shape noise in 
each survey. 

In addition, we will use simulations to estimate the uncertain¬ 
ties on our power spectrum measurements. Given a set of simula¬ 
tions containing both signal and noise components, we estimate the 
covariance matrix of the power spectrum estimates as: 

(AP b AP b ,) = <C P b ~ P b )(P v ~ P v )) mc, (34) 

where P b denotes the average of each band power over all simula¬ 
tions. 

Finally, we can also use simulations to validate our power 
spectrum analysis pipeline and to aid in the interpretation of null 
tests of the results. 


5.1 Signal simulations 

We generate the signal component of our simulations based on 
a ACDM cosmological model with parameter values taken from 
Planck Collaboration et al. (2014): Q m = 0.3175, cr 8 = 0.8347, 
H 0 = 67.1 km s -1 Mpc -1 , Cl b = 0.0486 and n s = 0.963, where 
Q m and Q b are the matter and baryon densities, cr 8 is the clustering 
amplitude in 8 h ~ l Mpc spheres, H 0 is the Hubble constant and n s is 
the spectral index of the primordial perturbations. We additionally 
assume a flat Universe, O a = 1 - Q m . 

To generate the input model spectra, we also need to model 
the redshift distributions of the two surveys. As mentioned earlier, 



Z 

Figure 8. Normalised galaxy redshift distributions adopted for the SDSS 
population (solid line) and for the FIRST population (dashed line). We use 
these distributions to generate the model shear power spectra for our simu¬ 
lations. 


using the S 3 simulation (Wilman et al. 2008) we have estimated 
the median redshift of the FIRST survey to be z RIRST = 1.2. For 
SDSS, we adopt a median redshift for the PHOTO catalogue of 
z sdss _ q 53 f rom Sypniewski (2014). This latter study utilised a 
method for estimating photometric redshifts using boosted decision 
trees and the known redshifts of the SDSS SPECTRO sample. We 
then model both the SDSS and FIRST redshift distributions using 
the parameterized form, 


n(z) = p 


exp 



(35) 


where p = 1.5, z* = z m /1.412 and z m is the assumed median red¬ 
shift (0.53 for SDSS and 1.2 for FIRST). The resulting redshift 
distributions are shown in Fig. 8 . 

Using the above parameter values and n(z) specifications, we 
generate the three model shear power spectra (the auto power spec¬ 
tra for the two surveys and the cross-power) according to equa¬ 
tion (21). For the 3-D matter power spectrum, we use the Bond & 
Efstathiou (1984) CDM transfer function. We also use the Halofit 
formula (Smith et al. 2003) to predict the non-linear P 5 (k , r), al¬ 
though for the scales of interest here, the non-linear contribution is 
negligible. 
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We use these power spectra to generate correlated Gaussian 
realisations of the shear fields in the two surveys. To create the 
correlated fields, at each multipole, £, we construct the 2 x 2 power 
spectrum matrix, C K f (l)K(J) where {/, j) denote the two surveys. Taking 
the Cholesky decomposition of this matrix at each multipole, L^ 7 , 
defined by 

Cf‘HA = '£ i L?Lf, (36) 

k 

we generate random Gaussian realisations of the spin-2 spherical 
harmonic coefficients of the two shear fields as 

a £0 ~ ^ Ljch, 

j 

(37) 

j 

where G l {m is an array of unit norm complex Gaussian random de¬ 
viates. The resulting fields transformed to real space via a spin- 
2 transform will then be correctly correlated according to the in¬ 
put auto and cross-power spectra. When generating the real-space 
shear fields, we use the same spatial resolution as was used in con¬ 
structing the shear maps from the real data (healpix A S j de = 64, 
corresponding to 0.85 deg 2 pixels). 

For each survey, we then create mock galaxy shape catalogues 
as follows. For each object in the real catalogues, we replace its 
measured ellipticity components with the components of the sim¬ 
ulated shear signal at the appropriate sky location in the healpix- 
generated maps. Note that we do not randomize the galaxy posi¬ 
tions. We choose to keep these fixed as we wish to mimic the real 
data as closely as possible in the simulations. One aspect of this is 
that we require the structure in the simulated galaxy number density 
maps to be the same as that for the real data (see Figs. 4 and 7). The 
observed large-scale structure in the number density maps is clearly 
dominated by observational selection effects rather than being due 
to intrinsic galaxy clustering. It is therefore appropriate to use the 
exact same number density fluctuations in the simulations in order 
to capture the impact of these selection effects. Using the unaltered 
galaxy positions from the real catalogues in all of the simulations 
ensures that this is the case. 


5.2 Noise simulations 

To create realisations of the random noise, we assign to each sim¬ 
ulated galaxy the measured ellipticity components of a real galaxy, 
sampled at random from the real galaxy catalogues. This procedure 
should model the random noise properties of the data (due to both 
the intrinsic dispersion in galaxy shapes and measurement errors 
in the shape estimation step) and facilitate the estimation of the 
noise bias term, (A^)mc- In principle, our procedure may slightly 
over-estimate the noise since the dispersion in the observed galaxy 
shapes will be enhanced slightly due to any lensing signal present. 
However, we expect this effect to have a negligible effect on our 
analysis - the ellipticity variance due to the combination of intrin¬ 
sic shape dispersion and measurement errors will greatly dominate 
over the lensing contribution for any plausible cosmological sig¬ 
nal. In any case, our cross-power spectrum measurements are not 
affected by the noise bias term, (A^)mc- 
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Figure 9. The y t component of the contamination template used to model 
the FIRST residual beam systematic discussed in Section 4.1.2, and dis¬ 
played in Fig. 2. The y r component of the contamination was assumed to 
be zero. 



5.3 Modelling small-scale systematic effects 

In our FIRST simulations, we also include a model of the spurious 
shear systematic, induced by the FIRST beam residuals, described 
in Section 4.1.2, and illustrated in Figs. 1 & 2. To model this effect, 
we first created a template of the observed spurious shear signal. 
This template is shown in Fig. 9 and has been normalised such 
that its azimuthally averaged tangential shear profile matches that 
measured from the real FIRST survey, shown in Fig. 1. 

We then use this template to model the contaminating influ¬ 
ence of each source detected in the FIRST survey on all other 
FIRST sources. The spurious signal induced on one source A, due 
to another source B, is modelled as 

e[ pur = y t (ARA, AS) cos(20), 

e pur = y t (ARA, AS) sin(2<9), (38) 

where yf(ARA, AS) is the tangential shear extracted from the tem¬ 
plate at the position (in {ARA, Ab}-space) of source A relative to 
source B. In equation (38), 6 describes the orientation of the great 
circle connecting the two sources with respect to the chosen refer¬ 
ence axis of the yi /2 coordinate system. 


5.4 Power spectra from simulations 

We create 100 signal-only simulated datasets as described in Sec¬ 
tion 5.1, each containing a known cosmic shear signal, appropri¬ 
ately correlated between the FIRST and SDSS realisations. Each 
realisation set also contains a noise-only simulation (Section 5.2) 
which is uncorrelated between datasets, and a model of the resid¬ 
ual FIRST beam systematic which is the same for each realisation 
(Section 5.3). For each realisation, the sum of all three of these 
components produces a full mock realisation of the FIRST and 
SDSS shear catalogues. 

We process each simulated catalogue in exactly the same way 
as we do for the real data. First, maps of the shear components 
are created from the simulated FIRST and SDSS ellipticities using 
simple averaging within A side = 64 healpix pixels. These maps 
are then multiplied by the appropriate weight maps (Figs. 4 & 
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Figure 10. The recovered mean SDSS-SDSS (left), SDSS-FIRST (centre) and FIRST-FIRST (right) Cf (blue squares), Cf (red circles) and cf (green 
diamonds) power spectra. The solid green curve shows the input Cf power spectra. The input cf and C'f power spectra were set to zero. Continuous and 
dashed line error bars show the error in the mean recovered values, and the uncertainty associated with a single realisation or measurement, respectively. 


7) before they are passed through the power spectrum estimation 
pipeline described in Section 2.2. 

Fig. 10 shows the mean recovered power spectra from the 100 
simulations. We estimate the three possible shear power spectra: 
the SDSS auto power, the FIRST auto-power and the cross-power 
between the two surveys. For each of these, in addition to the Cf 
(or £-mode) power spectra, we have also measured the Cf (or B- 
mode) power and the cross-correlation, cf (or EB). We measure 
the power spectra in four band powers (1<^<32;33<^<64; 
65 < i < 96; 97 < i < 128). For each band power we display 
two error bars: the larger (dashed) errors show the uncertainty as¬ 
sociated with a single realisation or measurement (calculated as 
crp b = (P 2 ) l/2 where (P 2 ) are the diagonal elements of the covari¬ 
ance matrix, equation 34) while the smaller errors show the error in 
the mean recovered band powers, o-p = crp b l a/A^mc- 

Inspecting the results, we see an unbiased recovery of the in¬ 
put Cf signal in all cases while the recovered Cf and Cf are con¬ 
sistent with zero, as required. The dashed error bars suggests that 
SDSS has the raw precision to detect the cosmic shear signal on 
these scales, the cross-power could potentially make a marginal 
detection while the FIRST dataset is too noisy on its own to de¬ 
tect the signal. Note that, on the relatively large scales of interest 
here, there is no noticeable effect from the residual beam system¬ 
atic that we have included in our FIRST simulations. However, we 
have confirmed with higher resolution simulations that this system¬ 
atic does bias the recovery of the FIRST auto-power spectrum, and 
increases the uncertainties on the FIRST-SDSS cross-power spec¬ 
trum, on smaller angular scales (at higher €). 


5.5 Recovery in the presence of large-scale systematics 

The results of the previous section demonstrate that our power 
spectrum analysis can recover a known input signal in the pres¬ 
ence of both random noise and the small-scale FIRST systematic. 
Those results also show that our analysis correctly accounts for the 
effects of the finite sky coverage and masking of the data. 

However, visual inspection of Figs. 3 & 6 suggests that there 
are additional large-scale systematics present in the data that are 
not modelled by our simulations. Our approach to dealing with 
these systematics is a central theme of this paper: by measuring 
the FIRST-SDSS cross-power spectrum, these systematics should 
drop out of the analysis, provided that the FIRST and SDSS sys¬ 
tematic effects are uncorrelated. To quantify the advantages of this 


approach, here we demonstrate the unbiased recovery of a simu¬ 
lated cross-power spectrum signal in the presence of either the real 
SDSS or the real FIRST large-scale systematic effects. 

To demonstrate unbiased recovery of the cross-power spec¬ 
trum in the presence of the SDSS systematics, we superimpose the 
SDSS signal-only simulations onto the real SDSS map. Along with 
any signals present in the real data (cosmic shear, noise and system¬ 
atics), the resulting coadded maps will now also contain a simulated 
signal that is correlated with that present in the FIRST simulations. 
We have passed these coadded “real + simulated” maps through 
our power spectrum estimation pipeline and the results are shown 
in Fig. 11. These recovered power spectra have undergone noise- 
bias subtraction and correction of the mask in exactly the same way 
as before. Fig. 11 shows that the recovery of the SDSS auto power 
spectrum is very heavily biased by the presence of unaccounted-for 
large-scale systematics in the data. 17 However, the right hand panel 
of Fig. 11 shows that the recovery of the SDSS-FIRST cross-power 
spectrum remains unbiased, despite the presence of the SDSS sys¬ 
tematics. Note that, although it is unbiased, the cross-power spec¬ 
trum is not completely unaffected by the presence of the SDSS sys¬ 
tematics: their presence results in a significant increase in the error 
bars of the cross-power spectrum recovery (compare the right-hand 
panel of Fig. 11 and the central panel of Fig. 10), due to chance 
correlations between the noise in the FIRST simulations and the 
SDSS systematic effects. We will account for these additional con¬ 
tributions to the errors in our final analysis. 

To assess the impact of the large-scale systematics in the 
FIRST shear maps, we have performed an equivalent analysis to 
what is described above but now we superimpose the FIRST signal- 
only simulations onto the real FIRST maps and cross-correlate 
these with the SDSS simulations. The results of this test are dis¬ 
played in Fig. 12 where we observe similar effects to what is de¬ 
scribed above - the FIRST auto-power spectrum is heavily biased 
by large-scale systematics while the cross-power spectrum remains 
unbiased. 


17 The alternative explanation - that the power observed in the left hand 
panel of Fig. 11 is of cosmological origin - is ruled out given current CMB 
and large-scale structure observations, and in any case is incompatible with 
the FIRST-SDSS cross-power spectrum that we measure later. The very 
large cf signal observed is a further strong indicator that the observed Cf 
signal in the SDSS auto power spectrum is due to large-scale systematics. 
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Figure 11. The recovered mean SDSS auto power spectra (left panel ) and SDSS-FIRST cross-power spectra (right panel ) from 100 simulations in the presence 
of large-scale systematic effects in the real SDSS shear maps. For a definition of the symbols see Fig. 10. Systematic effects in the SDSS data severely bias the 
recovery of the auto power spectrum, shown as the solid curve in the left hand panel. The cross-power spectrum recovery is unbiased in the presence of the 
same systematics, though the errors are amplified. 



Figure 12. The recovered mean FIRST auto power spectra (left panel) and SDSS-FIRST cross-power spectra {right panel) from 100 simulations in the presence 
of large-scale systematic effects in the real FIRST shear maps. For a definition of the symbols see Fig. 10. Systematic effects in the FIRST data severely bias 
the recovery of the auto power spectrum, shown as the solid curve in the left hand panel. As in Fig. 11, the systematics contribute additional errrors into the 
cross-power spectrum measurement though it remains unbiased. 


6 REAL DATA MEASUREMENTS 

In the previous section, we have demonstrated, using simulations, 
that we can extract an unbiased cosmic shear signal by cross- 
correlating the two datasets even in the presence of significant sys¬ 
tematic effects in both surveys, provided that the SDSS and FIRST 
systematics are not correlated with one another. We now apply the 
analysis to the real datasets. 

6.1 The SDSS-FIRST cross-power spectra 

To extract the cross-power spectrum from the real data, we apply 
the exact same procedure as was followed for measuring the power 
spectrum from the simulations in the previous section. Explicitly, 
we apply the power spectrum estimator of equation (26) to the real 
data maps, including both the noise debiasing step (based on the 
mean noise bias, (N Mc ), measured from noise only simulations) 
and the correction for the survey geometries and masks. 

To estimate errors, we once again use the scatter in the MC 
simulations (equation 34). However, in addition to the normal cos¬ 
mic variance and noise terms, we now also wish to include the en¬ 
hanced uncertainties due to chance correlations between random 


noise and large-scale systematics discussed in Section 5.5. To fa¬ 
cilitate this, we generate the total mock cross-power spectrum for 
the each simulation realisation as 

P l b - P b {sim l F x sim^) + P^sim^ x reals) + P^(real F x sim^), (39) 

where P b { sim^ x sim^) is the cross-power spectrum measured from 
the z th simulation set containing both signal and noise, P b (sim l F X 
reals) is the power spectrum measured by cross-correlating the / th 
FIRST simulation with the real SDSS map, and P b { real F x sim^) 
is the power spectrum measured by cross-correlating the z th SDSS 
simulation with the real FIRST map. The covariance matrix of our 
measurements can then be calculated from the scatter amongst the 
total mock power spectra of equation (39), and will now include un¬ 
certainties due to cosmic variance, random noise and the enhanced 
uncertainties due to the systematic effects. 

The cross-spectrum measured from the real FIRST and SDSS 
datasets is shown in Fig. 13. Comparing the real data measure¬ 
ment to the distribution of measured power spectra from simula¬ 
tions that include only noise and systematics, we find our measure¬ 
ment equates to a marginal detection of the C™ power spectrum 
at the 99% confidence level (a “~2.7cr detection”). The measured 
signal is in agreement (within ~lcr) with the model theory power 
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Figure 13. The cosmic shear cross-power spectrum measured from the 
SDSS and FIRST datasets. The blue squares show the measured E-mode 
signal {Cf), the red circles show the measured 5-mode {Cf) and the green 
diamonds show the measured EB correlation cf. The green solid curve is 
the expected Cf power spectrum for two populations with median redshifts 
of 0.53 and 1.2 respectively. 

spectrum where the median redshifts for the SDSS and FIRST pop¬ 
ulations are assumed to be 0.53 and 1.2 respectively. The 5-mode 
power spectrum (Cf) and the EB cross-correlation (Cf ) are both 
consistent with zero (within ~lcr). 

To further assess the significance of our measurement and the 
degree to which it is consistent with the expected signal in the 
concordance cosmological model, we have calculated a x 2 statis¬ 
tic from our suite of total simulated power spectra (equation 39) 
as: 

* 2 = “ p ?) 2 / <7 v < 4 °) 

b 

where P* is the expected value of the band power in the concor¬ 
dance model. We then compare the^ 2 value measured from the real 
data to the distribution of values from the simulated data sets. We 
have calculated this statistic for the three power spectra {Cf, Cf 
and Cf) separately and note that, for the latter two spectra, Pf = 0. 
The results of this test are shown in the top panels of Fig. 14. We 
see that our measurement from the real data is consistent with the 
simulation distribution, and hence with the expectation based on 
the concordance cosmology. 

The x 1 statistic tests for a particular type of discrepancy be¬ 
tween model and data. However, it does not test for all possible 
deviations. In particular, it it insensitive to the sign of P b - P*. We 
have therefore also performed a consistency test using the follow¬ 
ing statistic (which we loosely call^): 

X = E J (P b -Pf)/cr Pb . (41) 

b 

The results from this test are shown in the lower panels of Fig. 14 
and once again, reveal no inconsistencies between our measure¬ 
ments and the concordance cosmological model. 

6.2 Null tests 

We further assess the credibility of the data and results through a set 
of measurements designed to reveal any residual systematic s that, 
for whatever reason, might not be mitigated using our cross-power 
spectrum approach. For each of these null tests, we have processed 


each of our simulated datasets in exactly the same way in order 
to assign appropriate error bars to the null test measurements. The 
tests that we have performed are the following: 

(i) North-South. We split both the SDSS and FIRST data into 
North and South samples, with roughly the same number of sources 
in each sample. We then measure the power spectra for the two 
samples separately. Finally we subtract the signal between the two 
sets of spectra. 

(ii) East-West : This is the same as test (i) but splitting the data 
into East and West components. 

(iii) FIRST random : We randomly split the FIRST data into two 
subsets with the same number of galaxies. For each subset, a set of 
7 i and y 2 maps is then created. We take the difference of the shear 
maps for the two FIRST subsets and then estimate the cross-power 
spectrum of the resulting residual map with the SDSS shear map. 

(iv) SDSS random: The same as test (iii) but randomly splitting 
the SDSS data. 

(v) FIRST PSF: This test is designed to test whether the mea¬ 
sured signal is an artefact of the FIRST PSF. In addition to de¬ 
convolved shape estimates, the FIRST catalogue also includes es¬ 
timates of source shapes before PSF deconvolution. We replace 
our FIRST galaxy ellipticity catalogue with an ellipticity catalogue 
based on these uncorrected shape measurements for only those 
sources that are flagged as point sources (unresolved) in FIRST. 
The shear maps constructed from this catalogue should then trace 
the FIRST PSF accurately. We then estimate the cross-power spec¬ 
trum of the FIRST PSF map and SDSS shear catalogue. 

(vi) SDSS PSF: To test for artefacts associated with the SDSS 
PSF, we construct SDSS PSF maps using the quoted values in 
the SDSS survey for the reconstructed PSF at the location of each 
galaxy. The resulting PSF ellipticity maps are then cross-correlated 
with the FIRST shear catalogue. 

(vii) FIRST P(S) > 0.05: This final test is not strictly a null test. 
However, we implement it to test for any dependence of our mea¬ 
surement on the likelihood of the FIRST sources being sidelobe 
residuals. To perform this test, we create an alternate FIRST shape 
catalogue based on only those sources that have a higher chance of 
being a sidelobe, P{S) > 0.05. (Recall our main analysis is based 
on only sources with P{S) < 0.05.) We then repeat the FIRST- 
SDSS cross-power spectrum analysis using the P{S) > 0.5 alterna¬ 
tive FIRST catalogue. 

The Cf, Cf and Cf power spectra measured for all of these 
tests are shown in Fig. 15. In order to interpret these null tests, in 
Fig. Al, we once again plot the ^ 2 measurement from the real data 
alongside the x 2 histograms from the simulations. When calculat¬ 
ing the ^ 2 values for the null tests, we set Pf = 0 in equation (40). 
Examination of Figs. 15 and Al show that the data passes most of 
these null tests in the sense that the ^ 2 values from the real data are 
usually consistent with being a random sample of the simulation 
X 2 distribution. Potential problem cases are the East-West Cf and 
SDSS PSF Cf null tests. However, taken as a population, these tests 
do not reveal any major outstanding problems in the cross-power 
spectrum results. To further quantify the results of these tests, in 
Table 1, we list the probability to exceed (PTE) values for each test, 
which gives the probability that the ^ 2 value measured from the real 
data is consistent with being a random sampling of the simulation 
distribution. Low values of these PTE numbers indicate potential 
residual systematics. For completeness, in Fig. A2, we also present 
the measurements of ^ (equation 41 with Pf = 0) from the real 
data along with the simulation histograms. These tests reveal no 
obvious problems with the analysis. 
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Figure 14. Histograms showing the^ 2 values (equation 40; upper panels ) and^ values (equation 41; lower panels ) for the C™ (left), Cf (centre) and cf 
(right) power spectra as measured from the simulations. Over-plotted as the vertical blue line is the equivalent value for the real data measurements. 
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Figure 15. The C™ (left), cf (centre) and cf (right) power spectra measurements for the null tests described in Section 6.2. The displayed results are as 
follows: North-South : blue squares; East-West : red circles; FIRST random : green face up triangles; SDSS random : magenta horizontal lines; FIRST PSF: 
orange face left triangles; SDSS PSF : brown rectangles; FIRST P(S)>0.05: black crosses. 


Table 1. PTE values from the the^ 2 null tests described in Section 6.2. 


Null test 

r'KK 

w 

cf 

S~' K P 

North-South: 

0.64 

0.10 

0.10 

East-West: 

0.81 

0.01 

0.12 

FIRST random: 

0.41 

0.53 

0.52 

SDSS random: 

0.35 

0.10 

0.19 

FIRST PSF: 

0.96 

0.85 

0.87 

SDSS PSF: 

0.18 

0.67 

0.02 

FIRST P(S)>0.05: 

0.35 

0.31 

0.99 


6.3 Constraints on the FIRST and SDSS median redshifts, 
cr 8 and Q m 

Having demonstrated the validity of our results using null tests, we 
now use our cross-power spectrum measurement to constrain the 
median redshifts of the two populations, and the power spectrum 
normalisation, cr 8 and matter density, Q m . 

To constrain the median redshifts of FIRST and SDSS we per¬ 
form a grid-based likelihood analysis and generate the C K € K spec¬ 


trum for two populations with median redshifts in the range 0.05 < 
z s m DSS < 5.05 and 0.5 < z RIRST <5.5. The grid resolution in both 
directions is = 0.01. We fix the cosmological parameters val¬ 
ues at cr 8 = 0.8347 and Q m = 0.3175 (Planck Collaboration et al. 
2014). The model FIRST-SDSS cross-power spectrum is generated 
using equation ( 21 ) and is then averaged in band powers similarly 
to what was done during the power spectrum estimation. 

For each point in parameter space, we construct the x 2 misfit 
statistic (equation 40) which we then convert to likelihood values 
for a model with two degrees of freedom. The results are shown 
in Fig. 16. We immediately see that the constrains are symmetric 
about the z^ DSS = z RIRST line. However, for the purposes of this 
study we will make the reasonable assumption that z^ DSS < z RIRST . 
With this prior in place, the best fitting values for the median red¬ 
shifts are z^ DSS = 1.5 and z RIRST = 1.75. However, the data are 
also entirely consistent (within 1 cr) with the median redshifts de¬ 
rived from the literature of z^ DSS = 0.53 (Sypniewski 2014) and 
z RIRST = 1.2 (Wilman et al. 2008). 

We now fix the median redshifts for SDSS and FIRST to the 
literature values (zf n DSS =0.53 and z™ ST =1.2) and use our C K f mea¬ 
surements to place constraints on the cosmological parameters, cr 8 
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Figure 16. Joint constraints on the median redshifts of the SDSS and FIRST 
surveys obtained from fitting theoretical models to our C™ cross-power 
spectrum measurements. Cosmological parameters were kept fixed at the 
concordance values reported in Planck Collaboration et al. (2014). 



Figure 17. Joint constraints on the matter density, £l m and power spectrum 
normalisation, erg from fitting theoretical models to our C™ cross-power 
spectrum measurements. The median redshifts were fixed at z^ DSS = 0.53 
and zJJ, IRST = 1.2. 


and The results are shown in Fig. 17. The best fitting values 
for the two parameters are cr 8 = 1.5^ (68%) and Q m = 0.3 +q| 
(68 %). Our results are in agreement with the values in Planck Col¬ 
laboration et al. (2014) at the ler level. Our results are also in good 
agreement with the cr 8 constraint obtained from the FIRST cosmic 
shear analysis of Chang et al. (2004). 


7 CONCLUSIONS 

We have presented a cosmic shear cross-power spectrum analysis 
of the SDSS and FIRST surveys, a pair of optical and radio sky sur¬ 
veys with approximately 10,000 deg 2 of overlapping sky coverage. 

The motivation for our study has been to demonstrate the 
power of optical-radio cross-correlation analyses for mitigating 
systematic effects in cosmic shear analyses. The shear maps that 
we have constructed from both the SDSS and the FIRST cata¬ 
logues are severely affected by systematic effects. Measuring the 
auto-shear power spectrum from either of the datasets results in 
a heavily biased measurement - in both cases, the measured con¬ 
vergence power spectrum, C K [ is more than an order of magnitude 
larger than that expected in the concordance cosmological model. 
The presence of a large R-mode signal (Cf) in both the SDSS 
and FIRST auto power spectra is further evidence that the shear 
maps are contaminated by large-scale systematic effects. For our 
SDSS maps, these large-scale systematics dominate over the ran¬ 
dom noise on the scales of interest here ( £ < 120). For FIRST, the 
spurious power due to the large-scale systematics is approximately 
the same as the noise power spectrum. Note however that the ran¬ 
dom noise component in the FIRST shear maps is much larger than 
in the SDSS maps due to the much small galaxy number density in 
FIRST. Despite these very large shear systematics, we have demon¬ 
strated using simulations that one can still recover a cosmic shear 
signal using the FIRST-SDSS cross-power spectrum. Although the 
cross-power spectrum is not affected by uncorrelated systematics 
in the mean, the associated uncertainties on the cross-power spec¬ 
trum measurements are amplified by the presence of the FIRST 
and SDSS systematic effects. In our final analysis, our total errors 
on C k € k have increased by a factor of ~2.5 compared to those due to 


cosmic variance and random noise alone. Although not the focus 
of this paper, it is likely that this enhancement in the cross-power 
spectrum errors could be significantly reduced by applying more 
sophisticated techniques for correcting the galaxy shapes measure¬ 
ments for the effects of PSF anisotropy in both surveys, prior to 
construction of the final shear maps. 

By cross-correlating the SDSS and the FIRST data, we ten¬ 
tatively detect a signal in the C™ power spectrum that is incon¬ 
sistent with zero at the 99% confidence level. In contrast to al¬ 
most all cosmic shear analyses to date (which primarily probe the 
shear signal on sub-degree scales), our measurements constrain the 
weak lensing power spectrum on large angular scales. Our mea¬ 
surements probe the power spectrum in the fully linear regime, in 
the multipole range 10 < £ < 100, corresponding to angular scales 
2 ° < 6 < 20 °. 

Our results are consistent (within ~lcr) with the expected sig¬ 
nal in the concordance cosmological model, assuming median red¬ 
shifts of z s m DSS = 0.53 and z RIRST = 1.2 for the SDSS and FIRST 
surveys respectively. Our measurements of the odd-parity {Cf ) and 
parity-violating (Cf) cross-power spectra are fully consistent with 
zero, which demonstrates the success of the cross-power spectrum 
approach in mitigating systematic effects. We have also validated 
our analysis and results by performing a range of null tests on the 
data. 

We have used our CJf measurement to jointly constrain the 
median redshifts of the two surveys. With cosmological parameters 
fixed at their concordance values and assuming z^ DSS < z RIRST , we 
find best fitting value of z^ DSS = 1.5 and z RIRST = 1.75. However, 
our constraints are weak and our measurements are also consistent 
with the literature values of z^ DSS = 0.53 and z™ ST = 1.2. Fix¬ 
ing the median redshifts to these literature values, we have used 
our measurements to constrain the cosmological parameters Q m 
and cr 8 , where we find best fitting values of Q m = 0.30+q 2 and 
cr 8 = 1.5+gg (68% confidence levels), which are consistent with 
values quoted in Planck Collaboration et al. (2014) at the icr level. 

Although the detection of C K f presented here is tentative and 
lacks the precision acheived by state-of-the-art optical weak lens¬ 
ing surveys, we believe the analysis techniques developed in this 
paper will prove extremely useful for future high precision cosmic 
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shear analyses. In particular, we have successfully demonstrated 
that one can extract an unbiased cosmic shear signal even in the 
presence of severe shear systematics using the cross-power spec¬ 
trum approach. This type of analysis will be well suited for per¬ 
forming cross-correlation studies of future overlapping optical and 
radio surveys (e.g. with SKA, LSST and Euclid). Cosmic shear 
analyses of these future surveys will require exquisite control of 
systematic effects if they are to meet their science goals of per cent 
level constraints on dark energy and modified gravity thoeries. The 
cross-power spectrum approach that we have demonstrated for the 
first time in this paper represents a very promising tool for achieve- 
ing the required control of systematic effects. 
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APPENDIX A: RESULTS FROM NULL TESTS 

Fig. Al compares the^ 2 statistic (equation 40)) for the data with 
the distribution of x 2 measured from the simulations for the suite 
the null tests described in Section 6.2. Fig. A2 shows the same com¬ 
parison for the^f statistic of equation (41). 
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Figure Al. Histograms showing the distribution of x 2 values (equation 40) measured from the simulations for the suite of null tests described in Section 6.2. 
The results are shown for the C £ K (left column ), C £ (centre column) and C'f (right column) power spectra. Over-plotted as the vertiocal blue line is the 
equivalent value for the real data measurement. 
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Figure A2. Histograms showing the distribution of x values (equation 41) measured from the simulations for the suite of null tests described in Section 6.2. 
The results are shown for the C™ (left column ), (<centre column ) and C'f (right column ) power spectra. Over-plotted as the vertiocal blue line is the 

equivalent value for the real data measurement. 
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